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Abstract 

Background: The circadian clock enables living organisms to anticipate recurring daily and seasonal fluctuations in 
their growth habitats and synchronize their biology to the environmental cycle. The plant circadian clock consists of 
multiple transcription-translation feedback loops that are entrained by environmental signals, such as light and 
temperature. In recent years, alternative splicing emerges as an important molecular mechanism that modulates the 
clock function in plants. Several clock genes are known to undergo alternative splicing in response to changes in 
environmental conditions, suggesting that the clock function is intimately associated with environmental responses 
via the alternative splicing of the clock genes. However, the alternative splicing events of the clock genes have not 
been studied at the molecular level. 

Results: We systematically examined whether major clock genes undergo alternative splicing under various 
environmental conditions in Arabidopsis. We also investigated the fates of the RNA splice variants of the clock 
genes. It was found that the clock genes, including EARLY FLOWERING 3 (ELF3) and ZEITLUPE {ZTL) that have not 
been studied in terms of alternative splicing, undergo extensive alternative splicing through diverse modes of 
splicing events, such as intron retention, exon skipping, and selection of alternative 5' splice site. Their alternative 
splicing patterns were differentially influenced by changes in photoperiod, temperature extremes, and salt stress. 
Notably, the RNA splice variants of TIMING OF GAB EXPRESSION I {TOG!) and ELF3 were degraded through the 
nonsense-mediated decay (NMD) pathway, whereas those of other clock genes were insensitive to NMD. 

Conclusion: Taken together, our observations demonstrate that the major clock genes examined undergo 
extensive alternative splicing under various environmental conditions, suggesting that alternative splicing is a 
molecular scheme that underlies the linkage between the clock and environmental stress adaptation in plants. It is 
also envisioned that alternative splicing of the clock genes plays more complex roles than previously expected. 

Keywords: Arabidopsis thaliana, Circadian clock, Transcription factor, Alternative splicing, Nonsense-mediated decay 
(NMD), Environmental stress 



Background 

The circadian clock is an endogenous time-keeping system 
that coordinates the physiology and behavior of a living 
organism to its environment [1]. In plants, the clock 
modulates rhythmic leaf movement, elongation rate of 
hypocotyls, roots, and stems, stomata aperture, stem 
circumnutation, and flower opening [1,2]. 
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Three major interlocked feedback loops constitute the 
plant circadian clock: the central loop, the morning loop, 
and the evening loop [3-5]. The central loop is medi- 
ated by the reciprocal repression between the morning- 
phased MYB transcription factors, CIRCADIAN CLOCK 
ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPO- 
COTYL (LHY), and the evening-phased pseudo-response 
regulator TIMING OF CAB EXPRESSION 1 (TOC1) [6,7]. 
In the morning loop, CCA1 and LHY promote the tran- 
scription of PSEUDO-RESPONSE REGULATOR 9 (PRR9) 
and PRR7 genes [8,9]. Closing the loop, the PRR members 
inhibit the transcription of CCA1 and LHY genes by 
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sequentially binding to the gene promoters from early 
morning (PRR9) through mid-day (PRR7) to evening 
(PRR5) [10,11]. The evening loop is illustrated by TOC1 
and a hypothetical component Y, the expression of which 
is repressed by TOC1 and, in turn, activates TOC1 expres- 
sion [12]. Recent studies have shown that three evening- 
phased factors, EARLY FLOWERING 3 (ELF3), ELF4, and 
LUX ARRHYTHMO (LUX), form the EVENING COM- 
PLEX (EC), which represses PRR9 gene and LUX gene itself 
[13,14], indicating that the auto-inhibition of EC replaces 
the component Y in the evening loop [15]. 

The circadian system is substantially influenced by 
external cues. Phytochrome- and cryptochrome -mediated 
light signals mediate the induction of CCA1, LHY, and 
PRR9 genes [8,16,17]. Temperatures also affect the ampli- 
tudes and rhythms of the clock gene expression [18]. In 
addition, growth hormones and abiotic stresses modulate 
the clock function. It has been observed that accumulation 
of CC47, TOC1, and GIGANTEA (GI) gene transcripts is 
differentially regulated by abscisic acid, brassinosteroid, 
and auxin [19]. High light stress induces CCA1 gene [20], 
linking the clock with plant stress adaptation. 

The clock components are also regulated at the post- 
transcriptional and protein levels. It has been shown that 
the stability of CCA1 mRNA and the translation of LHY 
mRNA are influenced by light [21,22]. In addition, the 
F-box protein ZEITLUPE (ZTL) is responsible for the 
dark-induced degradation of TOC1 protein [23]. Fur- 
thermore, temperature-dependent phosphorylation of 
CCA1 modulates its binding to target gene promoters 
[24]. Most recently, chromatin remodeling and alterna- 
tive splicing of the clock genes have been described as 
fundamental processes in the regulation of the clock 
function [25]. 

Some of the clock genes have been shown to undergo 
alternative splicing in response to abiotic stresses in plants 
[26,27], among which temperature regulation of CCA1 
alternative splicing is best characterized. CCA1 alterna- 
tive splicing produces two protein isoforms, the full-size 
CCAla form and the truncated CCAlp form that lacks 
the MYB DNA-binding motif [27]. CCAlp competitively 
inhibits CCAla activity by forming nonfunctional hetero- 
dimers that are excluded from DNA binding. CCA1 alter- 
native splicing is suppressed by low temperatures. Under 
low temperature conditions, CCAlp production is reduced, 
and thus CCAla activity is elevated, leading to the stimu- 
lation of freezing tolerance [27], linking the clock with 
temperature response. 

Recently, it has been reported that alternatively spliced 
RNA isoforms of some clock genes are degraded through 
the nonsense-mediated decay (NMD) pathway [28-33], 
unlike the productive alternative splicing of CCA1 gene. 
NMD has evolved as an mRNA quality control mechan- 
ism that degrades mRNA molecules harboring premature 



termination codons (PTCs), which generate truncated 
proteins that are harmful to cellular energy metabolism, 
and those having aberrantly long 3' untranslated regions 
(3'-UTRs) [32,33]. It is thus possible that alternative 
splicing serves as a precise mechanism for controlling 
the mRNA levels of the clock genes, depending on en- 
dogenous and external conditions. 

In this study, we systematically investigated the alterna- 
tive splicing patterns of major clock genes under various 
environmental conditions. We also examined the fates of 
the RNA splice variants. Our study shows that alternative 
splicing of the clock genes is differentially influenced by 
photoperiod and a variety of abiotic stresses. The results 
of our study show that although RNA splice variants of 
some clock genes are predicted to encode truncated ver- 
sions of the authentic proteins, those of other clock genes 
do not appear to encode specific proteins and, instead, are 
degraded through the NMD pathway. It is envisioned that 
alternative splicing plays more complex roles in the clock 
function than previously expected. 

Results 

Major clock genes undergo extensive alternative splicing 

On the basis of the prevalence of alternative splicing 
events in the plant circadian clock genes in the literature 
[20,26,27,34,35], we anticipated that alternative splicing 
of the core clock genes constitutes a critical component 
of the clock function. Previous reports have shown that 
alternative splicing of CCA1 is suppressed by low tempera- 
tures [20,27,35]. The alternative protein isoform (CCAlp), 
which lacks the protein domain required for DNA binding, 
acts as a dominant negative regulator of the authentic 
CCA1 transcription factor (CCAla), thus providing a self- 
regulatory circuit that links the clock with temperature 
stress response. 

To extend our understanding of the functional relation- 
ship between the clock genes and environmental stress 
responses, we selected a group of major clock genes that 
constitutes the plant circadian clock and investigated 
whether these undergo alternative splicing and their alter- 
native splicing patterns are altered under environmental 
stress conditions. 

Analysis of gene structures deposited in the public 
databases and literature search revealed that PRR7, 
PRR9, TOC1, and ZTL genes as well as CCA1 gene 
undergo alternative splicing [26,27,34,35], each producing 
two or more RNA splice variants (Figure 1). For each 
clock gene, the a transcript represents the RNA splice 
variant that retains all the exons but do not have any in- 
trons. The transcript represents the one that exists at 
the highest level among the RNA splice variants other 
than the a transcript. 

CCA1 alternative splicing is mediated by the retention 
of intron 4 and introduces a PTC into CCA 1/3 transcript 
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Figure 1 Genomic structures of major clock genes. The clock 
gene sequences were analyzed using the softwares provided by the 
TAIR database. The predicted genome structures of CCA1 (A), PRR7 
(B), PRR9 (C), TOO (D), ELF3 (E), and ZTL (F) genes are displayed. 
Black boxes depict exons. White boxes are 5' and 3' UTRs. F and R 
are primers used for RT-PCR analysis of the RNA splice variants (see 
Figure 2). The a transcripts encode full-size, authentic proteins, and 
the jS transcripts encode truncated forms. Asterisks indicate premature 
termination codons (PTCs). b, bases. The genomic structure of CCA 1 
gene, which has already been reported by us [27], was included here 
for the benefit of the reader. 



(Figure 1A). PRR7 alternative splicing is somewhat com- 
plicated. It is mostly mediated by the retention of intron 3, 
resulting in PRR7/3 transcript (Figure IB and Additional 
file 1). A PTC is introduced into the PRR7/3 transcript. 
Notably, it is also mediated by the skipping of exon 4 and 
the retention of introns 2 and 3 [26,34,35]. PRR9 alter- 
native splicing is unique, among others, in that the 
major alternatively spliced variant (PRR9/3) is produced 
by selection of alternative 5' splice site in intron 2 



(Figure 1C and Additional file 2). The presence of two 
additional RNA splice variants has also been recently 
reported [26,34,35]. 

A single TOC1 cDNA sequence was identified in the 
TAIR database. However, it has been shown that an al- 
ternative splicing event occurs by the retention of intron 
4 [26,34], introducing a PTC into TOC1/3 transcript 
(Figure ID and Additional file 3). It has been reported 
that RNA splice variants of ELF3 gene are hardly detected 
in wild-type plants, but several RNA splice variants are 
detected in the skip-1 mutant, which is defective in its 
splicing machinery [34], possibly due to the retention 
of intron 2 or 3 (Figure IE). We found that the ELF3 
gene undergoes alternative splicing in wild-type plants 
(Additional file 4). In addition, it was found that the 
ELF3/3 transcript is derived from the inclusion of a new 
alternative exon and a PTC is introduced into the splice 
variant. 

There are two Zri-specific cDNA sequences (ZTLa 
and ZTLfi) in the public database. Sequence comparison 
and direct sequencing of RT-PCR products revealed that 
the ZTL alternative splicing is mediated by the retention 
of intron 2 (Figure IF and Additional file 5). The ZTLfi- 
encoded protein has been considered as an authentic 
ZTL enzyme in the literature [23], which is probably 
because the abundance of the ZTLft transcript is much 
higher than that of the ZTLa transcript (see below). 

The modes of splicing events are diverse in the clock genes 

The abundances of RNA splice variants other than a 
and ft transcripts were relatively very low in most cases 
([26,34,35], this study). We therefore decided to further 
investigate only the a and /? transcripts for each clock 
gene. The predicted alternative splicing modes of the clock 
genes were verified by cloning of the RNA splice variants 
by RT-PCR and direct DNA sequencing (Additional files 1, 
2, 3, 4, and 5). Total RNA samples were subjected to 
RT-PCR using primer pairs that are specific to each 
RNA splice variant. The results showed that all the RT- 
PCR products have the sizes that were inferred from the 
predicted alternative splicing modes of the clock genes 
(Figure 2A). No RT-PCR products were detected when 
reverse transcription was omitted prior to PCR amplifi- 
cations, indicating that total RNA samples used were 
not contaminated with genomic DNA. 

The modes of alternative splicing are diverse in the 
clock genes (Figure 2B). Retention of specific introns 
mediates the alternative splicing of CCA1, PRR7, TOC1, 
ZTL, and ELF3 genes. Exon skipping is involved in PRR7 
alternative splicing. Meanwhile, alternative 5' splice site 
contributes to PRR9 alternative splicing. Alternative spli- 
cing of ELF3 gene was the most complicated. Retention of 
intron 2 or 3 has been implicated in the ELF3 alternative 
splicing [34]. However, direct sequencing of PCR products 
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Figure 2 Detection of RNA splice variants of the clock genes. A. Detection of RNA splice variants by RT-PCR. Total RNA samples were isolated 
from 10-day-old Col-0 plants grown on MS-agar plates under LDs at peak ZT point for each clock gene and subject to RT-PCR. Gene-specific F 
and R primer sets, as indicated in Figure 1, were used to detect the transcript isoforms of each clock gene. PCR reactions were also performed 
without reverse transcription (-RT) to verify the lack of genomic DNA contamination. The sizes of the PCR products are provided at the bottom 
of the figure. SM, DNA size marker, bp, base pair. B. Modes of splicing events. Detection methods for the alternative splicing events are listed in 
the 3 rd column with the references indicated in the 4 th column. The nucleotide sequences of the RNA splice variants were determined (This work) 
or verified by direct DNA sequencing in this work. RNA-seq, RNA sequencing. Sanger, DNA sequencing by Sanger method. 



revealed that an additional RNA splice variant (ELF3/3), 
which is probably most abundant among the splice vari- 
ants, was produced by the inclusion of an alternative 
exon. 

We measured the absolute amounts of the RNA 
splice variants of each clock gene by qRT-PCR analysis 
(Figure 3A), as has been described previously [36,37]. 
Ten-day-old plants grown on MS-agar plates under 
long days (LDs, 16-h light and 8-h dark) were harvested at 
zeitgeber time (ZT) points of peak expression for individ- 
ual clock genes (e.g. ZTO for CCAl and ZTL, ZT8 for 
PRR9, ZT4 for PRR7, and ZT12 for TOC1 and ELF3), 
thereby maximizing the detection sensitivity of a small 
quantity of mRNA. Absolute quantitation of the a and (3 
RNA splice variants of each clock gene showed that the 
ratios (%) of p/a + (3 were variable among them (Figure 3B). 



The ratio of CCAl RNA splice variants was 34.32%, simi- 
lar to what has been previously reported [27]. Those of 
the RNA splice variants of PRR7 and PRR9 genes were ap- 
proximately 29%. In contrast, those of TOC1 and ELF3 
genes were relative low (<10%). One distinction was ZTL 
gene. The ft transcript level was higher than that of the a 
transcript (fila + ft = -75%), which was in contrast to the 
other clock genes. It is assumed that the physiological sig- 
nificance of alternative splicing varies in each clock gene. 

Some RNA splice variants of the clock genes are 
degraded by NMD 

Alternatively spliced RNA variants containing a PTC 
enter either the productive or unproductive pathway. In 
the productive pathway, the mRNA is translated into a 
protein that is structurally distinct from the authentic 
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Figure 3 Absolute quantification of alternatively spliced RNA 
variants. Ten-day-old Col-0 plants grown on MS-agar plates were 
used for the extraction of total RNA samples. To maximize the 
sensitivity of detection, plants were harvested at the phase of peak 
expression for each gene. A series of 10-fold dilutions of plasmid 
DNA containing each gene sequence was used to generate a standard 
curve. The regression line from the dilution curve was used to 
determine the concentration of each RNA splice variant. Black 
circles represent the absolute amounts of RNA splice variants 
(A). C T; threshold cycle. The percentages of /3/a + /3 were calculated for 
each clock gene (B). 



protein. One example is the alternative splicing of CCA1, 
in which the CCAlp protein isoform possesses protein 
domains required for dimer formation and transcriptional 
activation but lacks the MYB DNA-binding domain 
[27]. In contrast, in the unproductive pathway, the tran- 
script is degraded via the NMD- mediated degradation 
pathway [30]. 

To investigate the fates of the RNA splice variants of the 
clock genes, we employed two assay systems: cycloheximide 
(CHX) treatment and NMD-defective Arabidopsis mutants 
that are routinely employed for this purpose in the litera- 
ture. NMD requires translation, and thus the translational 
inhibitor CHX suppresses the NMD-mediated degradation 
of mRNA [38,39]. Wild-type Arabidopsis plants were 
treated with CHX, and the levels of ft transcripts were 
determined by qRT-PCR. It was found that whereas the 
levels of PRR7/3, PRR9/3, and ZTL/3 transcripts were not 
influenced by CHX treatments (Figure 4A, left panels), 
those of TOClfi and ELF3/3 transcripts were signifi- 
cantly elevated after CHX treatments (Figure 4B, left 
panels). 

We next examined the ft transcript levels of each clock 
gene in the upfl-5 and upf3-l Arabidopsis mutants, in 
which the NMD pathway is impaired [40,41]. The levels 
of PRR7/3, PRR9/3, and ZTL/3 transcripts in the mutants 
were comparable to those in wild-type plants (Figure 4A, 
right panels). In contrast, those of TOClfi and ELF3/3 
transcripts were higher by approximately two-fold in the 
mutants than in wild-type plants (Figure 4B, right panels). 
We also examined the levels of TOC1/3 and ELF3fi 
transcripts in the upfl-5 and upf3-l mutants under 
heat conditions. The ft transcript levels were even 
higher in the mutants than in wild type plants when 
grown at 37°C (Additional file 6). The more prominent 
differences in the TOClfi and ELF3fi transcript levels 
at 37°C is due to the heat-induced alternative splicing 
of TOC1 and ELF3 genes (see below). Based on these 
observations, it was concluded that whereas the PRR7fi, 
PRR9fi, and ZTLfi transcripts are likely to encode specific 
proteins, like the CCAlfi transcript [27], the TOClfi and 
ELF3/3 transcripts are probably targeted by NMD. The 
sensitivity of the TOC1/1 and ELF3fi transcripts to NMD 
is also consistent with the notion that the steady-state 
levels of NMD target mRNAs were relatively low in 
many cases [30,42,43]. 

Since the identities of ZTLa and ZTLfi transcripts are 
currently unclear ([23], this work), we also examined 
the effects of CHX and upfl-5 and upf3-l mutations on 
the accumulation of ZTLa transcript. We found that the 
ZTLa transcript level was not affected by CHX treat- 
ments (Additional file 7). It was also unaltered in the 
upfl-5 and upf3-l mutants, like that of ZTLfi transcript 
under identical assay conditions. It is therefore likely 
that the ZTLa transcript is not targeted by NMD 
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and, instead, encodes a distinct protein, like the ZTL/3 
transcript 



Protein isoforms of the clock components are defective in 
different functional domains 

Some RNA splice variants, such as PRR7ft, PRR9fi, and 
ZTIjS transcripts that are insensitive to NMD, were pre- 
dicted to encode truncated proteins that harbor altered 
protein structural organizations, as has been demon- 
strated with CCAlp protein isoform [27]. In many cases, 
these structural alterations in the truncated forms include 
deletions, insertions, or substitutions of certain protein 
domains [44]. 



We analyzed the structural organization of the predicted 
protein isoforms of CCA1, PRR7, PRR9, and ZTL using 
the analysis tools in the SMART and Pfam databases 
(http://smartembl-heidelberg.de/ and http://pfam.sanger. 
ac.uk/, respectively). The amino acid sequences of the 
protein isoforms were obtained either from the TAIR 
database or deduced from the nucleotide sequences of 
RT-PCR products. 

Two possible translation products were deduced from 
PRR7fi transcript. One protein isoform would be a trun- 
cated form containing the N-terminal pseudo-receiver 
(PR) domain, which is generated by the translation from 
the start codon to PTC (Figure 1). In this translation 
scheme, the PRR7/3 transcript harbors a long 3-UTR. It 
has been previously shown that alternatively spliced RNA 
variants having long 3'-UTRs are frequently targeted by 
NMD [30]. The other protein isoform is a truncated form 
lacking the N-terminal PR domain, which was marked 
as PRR7(3 (Figure 5A). On the basis of the structural 
similarity of PRR7(3 to CCAlp and PRR9|3 and the in- 
sensitivity of PRR7/3 transcript to CHX, we believe that 
the PRR7/3 transcript encodes the PRR7p protein that 
harbors the N-terminal truncation. 

PRR7p and PRR9p protein isoforms possess the CON- 
STANS, CONSTANS-LIKE, and TOC1 (CCT) domains 
but lack the N-terminal PR domain (Figure 5), which 
mediates interactions with other proteins, such as PRR5 
[23,45-47]. The overall structures of the predicted ZTLa 
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and ZTLp protein isoforms were similar to each other 
except for the short C-terminal sequence. The ZTLp 
isoform is slightly smaller than the ZTLa isoform by 17 
residues (Figure 5A). We were unable to identify any 
distinct protein motifs in the C-terminal region of the 
ZTL proteins, and thus it is currently unclear whether 
the two ZTL protein isoforms are functionally distinct 
or not. 

Our data showed that TOClfi and ELF3/3 transcripts 
were targeted by NMD and were not expected to encode 
any proteins (Figure 5B). 

Short days suppress the alternative splicing of T0C1 and 
ELF3 genes 

Plants use the circadian clock to monitor daylength 
changes in inducing seasonal developmental responses 
[48,49]. We therefore hypothesized that photoperiod 
influences the alternative splicing patterns of the clock 
genes. 

Arabidopsis plants were entrained to either LDs or 
short days (SDs, 8-h light and 16-h light), and the levels 
of alternatively spliced RNA variants were compared by 
qRT-PCR. The amplitudes and rhythms of CCA1 and 
PRR7 gene expression were not detectably altered under 
SDs (Figure 6A and B, respectively). Under SDs, PRR9 gene 
was induced, but its expression rhythms were maintained 



(Figure 6C). It seems that the alternative splicing of the 
morning-phased genes is not discernibly influenced by 
photoperiod. 

We observed that the levels of TOCla and ELF3a 
transcripts were higher under SDs than under LDs, evi- 
dently during the night (Figure 6D and E, respectively). 
In contrast, the levels of TOC1/3 transcript were lower 
during the night under SDs, and those of ELF3/3 transcript 
were not altered under identical conditions compared 
with LDs. Notably, the peak level of the TOClfi tran- 
script shifted from ZT12 under LDs to ZT8 under SDs. 
Together, these observations indicate that SDs sup- 
press the alternative splicing of the TOC1 and ELF3 
genes. There were no discernible effects of SDs on the 
alternative splicing of ZTL gene (Figure 6F). 

Low temperatures suppress CCA1 and ELF3 alternative 
splicing but induce TOC1 alternative splicing 

Low temperatures dampen the cyclic expression of the 
clock genes, resulting in the repression of the clock func- 
tion [18]. Similarly, low temperatures suppress the alterna- 
tive splicing of CCA1, and the imbalance between CCAla 
and CCAlp protein isoforms leads to disturbed circadian 
rhythms and induction of freezing tolerance [27]. It was 
therefore suspected that low temperatures would also 
affect the alternative splicing of other clock genes. 
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Figure 6 Effects of photoperiod on the alternative splicing of the clock genes. Ten-day-old Col-0 plants grown on MS-agar plates under 
either LDs or SDs were harvested at the indicated ZT points for the extraction of total RNA samples. The levels of the RNA splice variants of CCA1 
(A), PRR7 (B), PRR9 (C), TOO (D), ELF 3 (E), and ZTL (F) genes were determined by qRT-PCR. Biological triplicates were averaged. Bars indicate the 
standard error of the mean. 



Kwon et al. BMC Plant Biology 2014, 14:1 36 
http://www.biomedcentral.com/1471-2229/14/136 



Page 8 of 15 



Arabidopsis plants were exposed to 4°C, and the levels 
of alternatively spliced RNA variants were measured by 
qRT-PCR. To eliminate the effects of light- dark transi- 
tions, the assays were conducted under continuous light 
conditions. It was found that the rhythmic accumulation 
patterns of a and ft transcripts of the clock genes were 
significantly altered at low temperatures. CCA1 alterna- 
tive splicing was suppressed at low temperatures; the 
levels of CCA la transcripts were elevated, whereas 
those of CCA 1 ft transcripts remained low (Figure 7A), 
as previously described [27] . 

The levels of both PRR7a and PRR7/3 transcripts were 
lower during the subjective day and higher during the 
subjective night compared to those at 23°C (Figure 7B), 
indicating that low temperatures do not affect PRR7 
alternative splicing but diminish its rhythmic expression. 
PRR9 is functionally redundant with PRR7 [12]. However, 
the effects of low temperatures on PRR9 expression were 
distinct from those on PRR7 expression. The levels of both 
PRR9a and PRR9f$ transcripts were markedly elevated at 
low temperatures throughout the time course, and the 
rhythmic expression was enhanced (Figure 7C), indicating 
that low temperatures do not affect its alternative splicing. 

The effects of low temperatures on the alternative 
splicing of the evening-phased genes were quite diverse. 
While the levels of TOCla transcripts remained un- 
changed, those of TOC1/3 transcripts were markedly 



higher at low temperatures (Figure 7D), indicating that 
low temperatures induce TOC1 alternative splicing. 
The levels of ELF3a transcripts were largely unaffected 
but loosed rhythmicity at low temperatures (Figure 7E). In 
contrast, the levels of ELF3/3 transcripts were drastically 
reduced, showing that ELF3 alternative splicing is sup- 
pressed at low temperatures. ZTL expression was sup- 
pressed at low temperatures, but its alternative splicing 
remained unaltered (Figure 7F). 

Heat induces the alternative splicing of CCA1, PRR7, TOC1, 
and ELF3 genes 

Heat stress has become an important issue in the field 
because of recent global warming that extensively affects 
plant growth and distribution [50]. Because the clock is 
entrained at least in part by temperature, heat would 
certainly influence the clock function. However, little is 
known about the relationship between heat stress and 
the clock. We therefore examined the effects of heat on 
the alternative splicing of the clock genes. The heat assays 
were performed under continuous light conditions to 
eliminate the effects of light-dark transitions. 

Interestingly, the balance between a and ft transcripts 
varied among different clock genes, whereas most clock 
genes were induced at 37°C. The levels of CCA1/3 and 
PRR7/3 transcripts were significantly elevated at some ZT 
points after heat treatments (Figure 8A and B, respectively), 




Figure 7 Effects of low temperatures on the alternative splicing of the clock genes. Ten-day-old Col-0 plants grown on MS-agar plates 
under LDs were transferred to 4°C under continuous light conditions. Whole plant materials were harvested at the indicated ZT points. The levels 
of the RNA splice variants of CCA 7 (A), PRR7 (B), PRR9 (C), 70C7 (D), ELF3 (E), and ZTL (F) genes were determined by qRT-PCR. Biological triplicates 
were averaged. Bars indicate the standard error of the mean. 
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Figure 8 Effects of heat on the alternative splicing of the clock genes. Ten-day-old Col-0 plants grown on MS-agar plates under LDs were 
transferred to 37°C under continuous light conditions. Whole plant materials were harvested at the indicated ZT points. The levels of the RNA 
splice variants of CCA1 (A), PRR7 (B), PRR9 (C), TOO (D), ELF3 (E), and ZTL (F) genes were determined by qRT-PCR. Biological triplicates were 
averaged. Bars indicate the standard error of the mean. Note that the expression data at 23°C are identical to those in Figure 7. 



showing that their alternative splicing was accordingly 
induced. The levels of PRR9a and PRR9fi transcripts 
were elevated to a similar degree, showing that its alter- 
native splicing is not affected by heat (Figure 8C). The 
elevation of TOClft and ELF3/3 transcript levels were 
more prominent than that of TOCla and ELF3a tran- 
script levels (Figure 8D and E, respectively), suggesting 
that their alternative splicing was induced by heat. Heat 
effects were marginal on ZTL expression. The levels of 
both ZTLa and ZTL/3 transcripts were slightly elevated 
after heat treatments (Figure 8F). 

High salinity suppresses ELF3 alternative splicing 

Salt stress influences plant growth and developmental 
processes, such as flowering time, which is closely as- 
sociated with the clock function [51-53]. We therefore 
examined the effects of high salinity on the alternative 
splicing of the clock genes. 

It appeared that CCA1 and ZTL genes are not influ- 
enced by high salinity (Figure 9 A and F, respectively). 
Notably, PRR7 and TOC1 genes were suppressed by high 
salinity. The levels of both a and ft transcripts of these 
clock genes were reduced under high salt conditions 
(Figure 9B and D, respectively), showing that their alter- 
native splicing is not influenced by high salinity. ELF3 
gene was also suppressed by high salinity (Figure 9E), 
but the reduction of ELF3/3 transcript level was more 



prominent than that of ELF3a level, showing that ELF3 
alternative splicing is suppressed by high salinity. The 
levels of PRR9a and PRR9/3 transcripts increased with 
high salinity (Figure 9C), indicating that PRR9 gene is 
induced but its alternative splicing is not affected by 
high salinity. 

Discussion 

Effects of environmental conditions on the alternative 
splicing of the clock genes 

Rhythmic expression of stress response genes and distinct 
phenotypes of clock mutants under abiotic stress condi- 
tions underscore the close connection between the circa- 
dian clock and environmental stress response in plants. 
One of the best-understood mechanisms is the clock 
control of C-REPEAT BINDING FACTOR (CBF) genes 
that play a pivotal role in cold stress response [18,54]. The 
central oscillators CCA1 and LHY regulate the expression 
of the CBF genes by binding directly to their gene pro- 
moters [55]. The CBF genes are also directly regulated by 
PRR9, PRR7, and PRR5 [56,57]. In addition, the transcrip- 
tion of CBF target genes, such as COLD REGULATED IS 
A (CORISA) and RESPONSIVE TO DISSECATION 29 A 
(RD29A) [58,59], is clock-controlled [18]. 

Altered stress responses of various clock mutants further 
support the clock control of abiotic stress adaptation. The 
prr9 prr7 prrS triple mutants exhibit enhanced resistance 



Kwon et al. BMC Plant Biology 2014, 14:1 36 
http://www.biomedcentral.com/1471-2229/14/136 



Page 10 of 15 




ZTO 4 8 12 16 20 24 



ZTO 4 8 12 16 20 24 



ZTO 4 8 12 16 20 24 ZTO 4 8 12 16 20 24 



cl 30 



PRR9a 


-•-Mock 
-■-Salt 










ZTO 4 8 12 16 20 24 



ZTO 4 8 12 16 20 24 



32 1 



* 0 



TOCla -t- Mock 
[ -"-Salt 







ZTO 4 8 12 16 20 24 




ZTO 4 8 12 16 20 24 



ELF3a [^T 



- Mock 
-Salt 





ZTO 4 8 12 16 20 24 ZTO 4 8 12 16 20 24 



ZTLa 




ZTL/3 




ZTO 4 8 12 16 20 24 ZTO 4 8 12 16 20 24 



Figure 9 Effects of high salinity on the alternative splicing of the clock genes. Ten-day-old Col-0 plants grown on MS-agar plates under 
LDs were transferred to hydroponic MS medium containing 200 mM NaCI. The levels of the RNA splice variants of CCA 1 (A), PRR7 (B), PRR9 
(C), TOC1 (D), ELF 3 (E), and ZTL (F) genes were determined by qRT-PCR. Biological triplicates were averaged. Bars indicate the standard error 
of the mean. 



to drought and cold stresses [60]. TOC1 -deficient mutants 
display drought-tolerant phenotypes [61]. In addition, it 
has been shown that Arabidopsis plants that are defective 
in CCA1, LHY, and GI genes are susceptible to freezing 
temperatures [55,62]. 

Although the linkage between the clock and environ- 
mental stress responses has been widely explored, molecu- 
lar mechanisms and underlying signaling schemes have 
not been studied at the molecular level in most cases. It 
has been reported that low temperatures reduce the amp- 
litude of the clock gene expression [18]. Meanwhile, it is 
known that the clock genes are regulated by extensive 
alternative splicing, which is influenced by low tempera- 
tures. It is therefore evident that alternative splicing of the 
clock genes should be taken into the interpretation of the 
expression analysis data under abiotic stress conditions. 

This study shows that a group of major clock genes 
undergoes alternative splicing through a variety of spli- 
cing modes, such as intron retention, exon skipping, and 
selection of alternative 5 ' splice sites, resulting in two or 
more RNA splice variants for each clock gene. It was 
also found that photoperiod and abiotic stresses, such as 
temperature extremes and high salinity, broadly affect 
the alternative splicing of the clock genes. On the basis 
of the effects of CHX on the relative levels of RNA splice 
variants and expression assays in NMD -defective mutants, 



we propose that the alternative splicing of CCA1, PRR7, 
PRR9, and ZTL genes is productive with RNA splice 
variants encoding distinct proteins. In contrast, the 
RNA splice variants of TOC1 and ELF3 genes are predicted 
to be degraded through the NMD-mediated degradation 
pathway. 

Collectively, our data strongly support the notion that 
the major clock genes are also regulated at the posttran- 
scriptional level through alternative splicing in addition to 
the transcriptional control under both normal and envir- 
onmental stress conditions. Alternative splicing-mediated 
control of the clock genes would serve as a molecular 
scheme that incorporates environmental stress signals 
into the clock, as has been verified with CCA1 alterna- 
tive splicing that links low temperature signals with the 
clock [27]. 

In this work, we focused on two major RNA splice 
variants of each clock gene, although additional RNA 
splice variants have been identified or predicted for some 
of the clock genes examined (Figure 1). More works on 
the additional RNA splice variants are required to further 
extend our understanding on the linkage between alterna- 
tive splicing events of the clock genes and environmental 
stress responses. Searching for a full set of RNA splice 
variants of each clock gene, as has been performed by 
RNA sequencing method [20], will also be helpful for 
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the elucidation of the clock function in abiotic stress 
adaptation. We are currently working on plants that are 
impaired in the alternative splicing of each clock gene 
and those expressing a specific RNA splice variant to in- 
vestigate the physiological roles of the alternative splicing of 
the clock genes. 

Function of alternative protein isoforms 

Recent studies have shown that protein isoforms that lack 
specific functional domains, which are produced by the 
alternative splicing of transcription factor genes, act as 
competitive inhibitors of the authentic transcription factors 
by forming nonfunctional heterodimers [27,63]. The best- 
characterized mechanism is the dominant negative regula- 
tion of the CCA1 transcription factor (CCAla) by the 
protein isoform CCAlp. While the CCAlp isoform pos- 
sesses protein domains required for dimmer formation and 
transcriptional activation, it lacks the MYB domain neces- 
sary for DNA binding [27]. Therefore, CCAlp is capable of 
interacting with CCAla, forming CCAla-CCAlp heterodi- 
mers that are excluded from DNA binding. 

According to the domain structures of the protein iso- 
forms encoded by the NMD-insensitive RNA splice vari- 
ants, the PRR7p and PRR9p protein isoforms are predicted 
to function in a way that is distinct from that of CCAlp. 
Unlike CCAlp that lacks the MYB DNA-binding domain, 
PRR7p and PRR9p have the CCT domain, which is respon- 
sible for DNA binding, but lack the PR domain that medi- 
ates protein-protein interactions [7,23,45-47]. A plausible 
working mechanism of the PRR7p and PRR9p protein iso- 
forms would be that they compete with the authentic 
PRR7a and PRR9a transcription factors for binding to the 
promoters of target genes, as has been previously proposed 
[64]. Further investigations are required to determine the 
functional modes of PRR7p and PRR9p. 

Two RNA splice variants of ZTL gene are also insensi- 
tive to NMD, and two protein isoforms, ZTLa and ZTLp, 
are expected to be produced. The ZTLa and ZTLp iso- 
forms are identical except for the far C-terminal se- 
quences; the former is larger than the latter by 17 
residues. The functional mode of ZTLp thus might differ 
from the p protein isoforms of other clock components. 
The lack of the C-terminal extension might also influence 
the protein conformation of the ZTLp isoform, which 
would affect its substrate specificity or enzymatic activity. 
The smaller ZTLp isoform has been annotated as the 
authentic ZTL protein in the literature [23,65]. We found 
that the level of ZTLfi transcript is much higher than that 
of ZTLa transcript, which is in contrast to the alft ratios 
of other clock genes. It is currently unclear whether ZTLa 
or ZTLp or both is an authentic enzyme. Phenotypic and 
physiological assays on transgenic plants that specifically 
express either ZTLa or ZTLfi cDNA would help clarify 
this uncertainty. 



NMD-mediated control of the clock gene expression 

Unlike the NMD-insensitive /? transcripts of CCAi, PRR7, 
PRR9, and ZTL genes, TOC1/3 and ELF3/3 transcripts 
are apparently targeted by NMD. The TOC1/3 and ELF3/1 
transcripts possess sequence features that are frequently 
observed in NMD substrates, in which they have 
splice junctions downstream of the PTC and very 
long 3'-UTRs [30]. 

Physiological roles of the NMD pathway are somewhat 
controversial. According to the "noise" hypothesis, NMD- 
sensitive RNA splice variants occur as a result of splicing 
error and are eventually removed through the NMD 
pathway [66]. In contrast, in the "regulated unproductive 
splicing and translation (RUST)" hypothesis, alternative 
splicing is coupled with NMD as a regulatory mechanism 
for monitoring the abundance of full-size RNA splice vari- 
ants [67]. We believe that the RUST hypothesis fits well 
into the alternative splicing of the TOC1 and ELF3 genes, 
based on the following reasons. First, the RUST hy- 
pothesis depicts that alternative splicing occurs through 
distinct modes of splicing events [26,34]. We found that 
the alternative splicing of the TOC1 and ELF3 genes is 
mediated by the retention of specific introns, supporting 
the notion that their alternative splicing is a regulated 
process rather than a simple splicing error. Second, their 
alternative splicing is regulated by environmental factors 
in a discrete manner. Production of the ELF3/3 transcript 
is suppressed by cold and high salinity conditions but 
induced under heat stress conditions. Third, whereas their 
alternative splicing is markedly influenced by abiotic 
stresses, the levels of TOCla and ELF3a transcripts are 
less affected under identical conditions. However, it is 
still possible that the RNA splice variants may be at 
least in part translated into proteins. It has recently been 
reported that some NMD targets are stabilized and trans- 
lated into proteins under certain conditions [68] . 

Conclusions 

We investigated the alternative splicing events of major 
clock genes under various environmental conditions and 
the sensitivity of their RNA splice variants to NMD. Alter- 
native splicing patterns of the clock genes were differently 
affected by changes in photoperiod and abiotic stresses, 
such as cold, heat, and high salinity. Based on the results 
of this study, we propose that alternative splicing of the 
clock genes, either by producing truncated isoforms 
that act as self-regulators or by regulating the abun- 
dance of full-size transcripts at the posttranscriptional 
level, contributes to the precise regulation of the clock 
function, particularly under fluctuating environmental 
conditions. It may also serve as a web that integrates 
environmental stress signals into the clock, providing 
an adaptation strategy in response to unpredictable en- 
vironmental changes. 
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Methods 

Bioinformatics software 

Gene sequences and their exon-intron structures were 
obtained from the Arabidopsis Information Resource 
(TAIR, http://www.arabidopsis.org/). Alternative splicing 
modes of PRR7 and T0C1 genes, which have not been an- 
notated in TAIR, were predicted based on the sequence 
analysis and the previous reports describing the predicted 
types of alternative splicing [26,34]. For ELF3 gene that has 
not been studied, the presence of alternatively spliced RNA 
variants was verified by direct sequencing of RT-PCR 
products. Protein domain structures were predicted using 
the SMART (http://smart.embl-heidelberg.de/) and Pfam 
(http://www.sanger.ac.uk/Software/Pfam/) databases. 

Plant materials and growth conditions 

Arabidopsis thaliana ecotype Columbia-0 (Col-O) was 
used, unless otherwise specified. The upfl-5 and upf3-l 
mutants, which have been previously described [26,34], 
were kindly provided by Dr. Jeong Sheop Shin (Korea 
University, Seoul, Korea) and Dr. Hee-Jeong Jeong (Kyung 
Hee University, Yongin, Korea). Plants were grown on X A 
X Murashige & Skoog media containing 0.6% (w/v) agar 
(hereafter referred to as MS -agar plates) in a growth 
chamber set at 23°C with relative humidity of 60% under 
either long day conditions (LDs, 16-h light and 8-h dark) 
or short day conditions (SDs, 8-h light and 16-h dark) 
with white light illumination (120 \iM photons m" 2 s" 1 ) 
provided by fluorescent FLR40D/A tubes (Osram, Seoul, 
Korea). 

Analysis of gene transcript levels 

Extraction of total RNA samples from appropriate plant 
materials and RT-PCR conditions have been described 
previously [69]. Total RNA samples were extensively 
pretreated with an RNase-free DNase to eliminate con- 
taminating genomic DNA prior to analysis. 

Quantitative real-time RT-PCR (qRT-PCR) was employed 
to determine the levels of gene transcripts. RNA sam- 
ple preparation, reverse transcription, and quantitative 
PCR were conducted according to the rules that have 
been proposed to ensure reproducible and accurate 
measurements [70]. 

qRT-PCR reactions were performed in 96-well blocks 
using an Applied Biosystems 7500 Real-Time PCR System 
(Foster City, CA) using the SYBR Green I master mix in a 
volume of 20 The PCR primers were designed using 
the Primer Express Software installed in the system and 
listed in Additional file 8. The two-step thermal cycling 
profile used was 15 s at 94°C and 1 min at 68°C. The 
eIF4A gene (At3gl3920) was included in the reactions 
as internal control to normalize the variations in the 
amounts of cDNA used [71]. All qRT-PCR reactions 
were performed in biological triplicates using RNA samples 



extracted from three independent plant materials grown 
under identical conditions. The comparative AAC T method 
was used to evaluate the relative quantities of each ampli- 
fied product in the samples. The threshold cycle (C T ) was 
automatically determined for each reaction using the 
default parameters of the system. The specificity of PCR 
reactions was determined by melt curve analysis of the 
amplified products using the standard methods installed 
in the system. 

Absolute quantification of gene transcripts 

Absolute quantification of gene transcripts was conducted 
as previously described [27]. The cDNAs of alternatively 
spliced RNA variants were subcloned into a pGADT7 vec- 
tor (Clontech, Mountain View, CA), and the absolute 
standard curve of each transcript was obtained by a series 
of 10-fold dilutions covering from 10" 19 to 10" 23 mol/ul, 
as described elsewhere [36,37]. Quantitative RT-PCR was 
conducted using a SYBR Green I master mix (Applied 
Biosystems) with splice variant- specific primers listed in 
Additional file 8. 

Abiotic stress treatments 

Arabidopsis plants grown for 10 days on MS-agar plates 
under LDs were used for abiotic stress treatments. For 
cold and heat treatments, plants were incubated at 4°C or 
at 37°C under continuous light conditions for appropriate 
time durations before harvesting plant materials. To exam- 
ine the effects of high salinity on the alternative splicing 
of the clock genes, plants were transferred to MS liquid 
medium containing 200 mM NaCl under continuous 
light conditions for appropriate time durations. 

Cycloheximide (CHX) treatments 

The CHX treatments were performed as described else- 
where [29,30]. Ten-day-old plants grown on MS-agar 
plates were transferred to MS liquid medium containing 
20 (iM CHX. Following vacuum infiltration for 10 min, 
the plants were incubated for 5 h at 23°C under normal 
growth conditions before harvesting plant materials for 
total RNA extraction. 

Additional files 



Additional file 1: Nucleotide sequence comparison of PRR7 gDNA 
and PRR7Q cDNA. The nucleotide sequence of PRR7f3 cDNA was 
determined by DNA sequencing of RT-PCR product and aligned with 
PRR7 genomic DNA {PRR7 gDNA) using the ClustalW software (http:// 
www.ebi.ac.uk/tools/msa/clustalw2/). Part of the aligned sequences 
containing exons 1, 2, 3, and 4 and introns 2, 3, and 4 was displayed. 
Intron 3, which is retained in the PRR7(3 transcript as a result of alternative 
splicing, is underlined (blue). A PTC (premature termination codon) is 
introduced into the PRR7(3 transcript (red asterisk). 

Additional file 2: Nucleotide sequence comparison of PRR9 gDNA 
and PRR90 cDNA. The nucleotide sequence of PRR9f3 cDNA was 
determined by DNA sequencing of RT-PCR product and aligned with 
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PRR9 gDNA using the ClustalW software. Part of the aligned sequences 
containing exons 1, 2, 3, 4, and 5 and introns 1, 2, 3, and 4 was displayed. 
The alternative splice site, which is used to produce the PRR9(3 transcript, 
is underlined. Sequence analysis revealed that the PRR9(3 transcript occurs 
by the alternative 5' splice site in intron 2 (blue). 

Additional file 3: Nucleotide sequence comparison of T0C1 gDNA 
and TOCIfi cDNA. The nucleotide sequence of TOC7/3 cDNA was 
determined by DNA sequencing of RT-PCR product and aligned with 
JOG gDNA using the ClustalW software. Part of the aligned sequences 
containing exons 4, 5, and 6 and introns 3, 4, and 5 was displayed. The 
retained intron 4, which is included in the TOC7/3 transcript as a result of 
alternative splicing, is underlined (blue). A PTC is introduced into the 
TOC7/3 transcript (red asterisk). 

Additional file 4: Nucleotide sequence comparison of ELF3 gDNA 
and ELF3& cDNA. The nucleotide sequence of ELF3(3 cDNA was 
determined by DNA sequencing of RT-PCR product and aligned with 
ELF3 gDNA using the ClustalW software. Part of the aligned sequences 
containing exons 2 and 3 and intron 2 was displayed. The alternative 
exon within intron 2, which is included in the ELF3(3 transcript as a result 
of alternative splicing, is underlined. Red boxes indicate conserved 'GT' 
and 'AG' sequences at the 5' and 3' ends of introns. Sequence analysis 
revealed that the ELF3(3 transcript occurs by the inclusion of an alternative 
exon consisting of 178 nucleotides within intron 2 (blue). A PTC is 
introduced into the ELF3f3 transcript (red asterisk). 

Additional file 5: Nucleotide sequence comparison of ZTL gDNA 

and ZTLa and ZTLfi cDNAs. The nucleotide sequences of ZTLa and ZTLfi 
cDNAs were determined by DNA sequencing of RT-PCR products and 
aligned with ZTL gDNA using the ClustalW software. Part of the aligned 
sequences containing exons 2 and 3 and intron 2 was displayed. Intron 
2, which is retained in the ZTLfi transcript as a result of alternative 
splicing, is underlined (blue). A PTC is introduced into the ZFLf3 transcript 
(red asterisk). The 3' untranslated region of the Z7L/3 transcript is shown 
in gray. 

Additional file 6: The fate of TOCIfi and ELF30 transcripts under 
heat stress conditions. Ten-day-old Col-0 plants and up/7 -5 and upf3- 7 
mutants grown on Vi X Murashige & Skoog media containing 0.6% (w/v) 
agar plates (hereafter referred to as MS-agar plates) were transferred to 
37°C for 12 h before harvesting whole plant materials for the extraction 
of total RNA. Levels of TOC7/3 and ELF 3$ transcripts were determined by 
quantitative real-time RT-PCR (qRT-PCR). Biological triplicates were averaged 
and statistically treated using Student f-test (*P<0.01). Bars indicate standard 
error of the mean. 

Additional file 7: The fate of ZTLa transcript. Plants were grown on 
MS-agar plates for 10 days under normal growth conditions. Ten-day-old 
Col-0 plants were transferred to liquid MS culture containing 20 uM 
cycloheximide (CHX). Following vacuum infiltration for 10 min, the plants 
were incubated for 5 h at 23°C under normal growth conditions before 
harvesting whole plant materials for the extraction of total RNA (left 
panel). The upfl-5 and up/3- 7 mutants were not treated with CHX (right 
panel). Levels of ZFLa transcript were determined by qRT-PCR Biological 
triplicates were averaged. Bars indicate standard error of the mean. 

Additional file 8: Primers used in qRT-PCR, RT-PCR, and gene cloning. 

F, forward primer; R, reverse primer. 
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